Hotspot Emission from a Freely Precessing Neutron Star 
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ABSTRACT 



Recent observations of IE 161348-5055 the neutron-star candidate at the center 
CN ■ 

J> ■ of the supernova remnant RCW 103, show that a component of its emission varies 

o 

00 ... 
t^J- , what one would expect for a freely precessing neutron star with a spin period of about 

one second. We produce light curves for a freely precessing neutron star with a hotspot. 

| By a suitable choice of parameters, we obtain light curves which are constant with 

rotational phase when the flux from the star reaches a maximum. At other phases of 

the precession, the flux varies as the star rotates but the total flux decreases by a factor 

of several. These models can explain the behavior observed from IE 161348-5055 and 

predict that the spin period should be detectable at minimum flux from sufficiently 



1. Introduction 



c/2 . sensitive measurements. 
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Garmire et al. (2000) have noted that the x-ray source, IE 161348-5055, located near the center 
of the supernova remnant RCW 103, has a sinusoidal light curve with a period of approximately 
six hours. The authors suggest that this period implies that the x-ray source has a low-mass 
companion with a six hour orbital period. In this paper, we examine an alternative possibility, that 
IE 161348-5055 is a freely precessing neutron star with a precession period of about six hours. 

The varying spectrum from IE 161348-5055 is similar to a blackbody yielding an effective 
area of a small fraction of a square kilometer (Garmire, private communication), much less than 
the surface area of a neutron star; therefore, the situation could be well approximated by a point 
source (a hotspot) on the surface of a freely-precessing neutron star. 

Free precession has often been invoked to explain long-period variability in neutron stars 
and neutron-star systems. Brecher (1975) attributed the 35-day cycle of Her X-l ( [Tananbaum 



et al. 1972p to the free precession of the neutron star secondary. |Jones (1988)| explained timing 
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residuals in the Crab pulsar as arising from small amplitude free precession of the neutron star. 
Cadez & Galicic (1996) found a small amplitude modulation in the optical flux from the Crab 



pulsar which Cadez, Galicic & Calvani (1997) cite as evidence of free precession. 

In §|2] we outline the kinematics of free precession, and in §||, we review the equations that 
determine the trajectory of light leaving the surface of the neutron star. §|] presents the light curves 
as a function of precessional and rotational phase for Newtonian, relativistic and ultracompact 
neutron stars with a emission from a hotspot. Finally, ^ places the results in a greater context. 



2. Free Precession 

A body will precess freely if its rotation axis does not coincide with one of its principal axes. 
More specifically, if the body is only slightly prolate or oblate, the angular velocity vector (c3) of 
the star will make a constant angle (k) with one of the principal axes (the 3— axis), forming the 
body cone, and will trace a cone in space (the space cone) with half-opening angle k as well. The 
rate of the precession is given by (|Goldstein 198CQ 



= Js T Tl uj 3 = eu 3 . (1) 
h 

Values of e = 10~ 3 — 10 -4 agree with the glitching behavior of neutron stars and the inferred 
gravitational-radiation spindown from the Crab pulsar ( jBhapiro Teukolsky 1983| ). Forced 
precession by an orbiting companion typically has a frequency lower by a factor of 2/(3Sl*u;3) 
where f2* is the angular frequency of the orbit. 



Internal magnetic fields may distort a neutron star significantly. Ostriker fc Gunn (1969) 
estimate the distortion of a neutron star due to internal fields, 

e^4xl0- 6 (3( J B p 2 15 )-( J B2 15 )) (2) 

where -B15 is the value of the magnetic field in units of 10 15 G and (...) denotes a volume- weighted 
average over the star. Here, B p and refer to the poloidal and azimuthal components of the 
magnetic field, respectively. One expects the internal fields of the neutron star to be significantly 
larger than the field inferred by magnetic dipole radiation due to the contribution of higher 



multipoles and the concentration of magnetic flux in flux tubes (e.g. Pines fc Alpar 1985 ) 



We would like to calculate the light curve from a hotspot on the surface of the freely precessing 
neutron star. Let us take the center of the cone that the angular velocity moves along to be the 
z— axis. Our line of sight (O) makes an angle £ with this axis and forms a plane with z— axis. We 
measure the phase of the precession (4>) relative to intersection of the space cone with this plane. 

The plane containing z and O also intersects the body cone. The angular momentum of the 
star points along the z-axis. In the body frame, the angular velocity of the star traces a cone 
centered on a principal axis of the star. We call this principal axis 3. The hotspot is located at fl 
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which makes an angle (5 with 3. Let us freeze the precession and the rotation of the star when 3 
points along z (see Figure [j]). At this orientation, the angle between the 3 — O plane and the 3 — // 
plane is 7. 

Using spherical trigonometry (see Figure |l|) yields the angle between the line of sight and the 
rotation axis, (((ft), an d the angle between the hotspot and the rotation axis a(<ft). 

cos (((ft) = cos k cos £ -|- sin n sin ( cos (ft (3) 
cosa((ft) = cos ft cos f3 + sin fcsin/3cos(0 — 7) (4) 



The star rotates as well as precesses. The phase of the rotation is given by rj. When rj = the 
hotspot lies in the u — O plane. The angle between the line of sight and the hotspot is 9 and is 
given by 

cos 9 = cos a((ft) cos (((ft) + sin a((ft) sin (((ft) cos 77 (5) 

If the hotspot emits isotropically, and we neglect gravitational lensing, the observed flux from the 
hotspot is simply proportional to cos# for \9\ < ir/2 and zero otherwise. 



3. Gravitational Lensing 

I □ 

Page (1995) presents a detailed treatment of the gravitational lensing of the surface of a neutron 
star. Since the light trajectory is bent, the zenith angle of our detector (5) as seen from the hotspot 
is no longer equal to 9. They are related by 

9 ^ = Jo TP - 2y)y - (1 - 2u)u 2 x' 2 (6) 

where y = GM/Rc 2 and x = sinS. For y < 1/3, the image of the hotspot i will be visible if 
9i + 2irj < 9(1). The flux from the hotspot is proportional to 



x(9i + 2irj) dx 
y sm(9i + 2vrj) 



(7) 

e=9i+2TTj 



such that \9i + 2vrj| < 9(1). 



For y > 1/3 the surface of the neutron star lies below the circular photon orbit; therefore, each 
hotspot yields an infinite number of images ( [Shapiro fc Teukolsky 19831) : 

lim 9(x;y > 1/3) = 00. (8) 

where 

£max = 3\/3y V 1 - 2 V- (9) 
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Fig. 1. — Definitions of various angies 
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4. Light Curves 

To construct a light curve, we must specify several angles (£, k, (3 and 7). During some portion 
of the precession period, the observed flux will be constant with orbital phase, if either f3 = k or 
£ = k. In the first case, the angular velocity vector will coincide with the location of the hotspot 
on the star. In the second case, the angular velocity vector will point along the line of sight once 
during each precession. 

4.1. (3 = K = m° 

To maximize the flux during the portion of the precessional period where the flux does not 
vary with rotational phase, we take 7 = and /3 = k = 90°, and to minimize the flux during the 
rest of the precession, we take £ = #(1). If 7 7^ 0, the maximum flux will not occur during that 
portion of precession when the flux does not vary with rotational phase. Taking (3 = k 7^ 90° will 
reduce the maximum flux, and £ 7^ #(1) will change the portion of the time when the hotspot is 
not visible. 

Figure ^ shows the observed flux from the hotspot as a function of the star's precessional 
and rotational phase. Only half of a precessional period is depicted. The flux varies at twice 
the precession rate. As Figure |3] shows, the mean flux over the rotational period varies nearly 
sinusoidally. Twice during the precessional period, when the mean flux reaches its maximum, the 
flux does not vary with the rotational phase of the star. At other stages of the star's precession, 
the hotspot spends much of the rotational period hidden behind the horizon on the neutron star 
surface. If (3 = n 7^ 90°, one finds that the mean flux varies at the precessional frequency, resulting 
in a light curve similar to that presented in Figure ||. 

4.2. £ = k = 9(1) 

In this case, we attempt to maximize the flux during the portion of the precession when the 
flux does not vary over the rotation. To do this, we take /? = £ = «; = 9(1) and 7 = 0. If the bending 
of the photon trajectories is neglected (i.e. as y approaches zero), 9(1) approaches 90° and this case 
reduces to the previous one. However, for a realistic neutron star with y ~ 0.2, we have new light 
curve which varies at the precession rate (not twice that rate as the previous case). Taking f3 7^ k 
or 7 7^ also yields a portion of the light curve when the flux does not vary with rotational phase, 
but this does not coincide with the period when the mean flux reachs its maximum. 
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Fig. 2. — Light curves as a function of precessional phase (</>) and rotational phase (rj) for y = 0, 0.2 
or R = oo and 11 km for a 1.4 M© neutron star. For both models, k = (5 = 90°, 7 = and £ = 0(1). 
A hotspot on a Newtonian star directly below the observer has a flux of unity in these units. 




Fig. 3. — Mean flux as a function of precessional phase (4>) for y = 0, 0.2 or R = oo and 11 km for 
a 1.4 Mq neutron star. The dashed curve gives the best-fitting function of the form a cos 2(f) + b. 
The geometry is as above. 
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4.3. Random Geometry 

The choices of £, 7, /5 and k that we have made previously are not generic. One would expect 
the values of £ and 7 from a particular neutron star to be random. The values of (5 and k are 
intrinsic to the star; therefore, one may find a physical motivation for a particular distribution of 
their values. Figure |H| presents two light curves for two randomly selected geometries. 

We created a random sample of several hundred geometries and found that approximately 
three percent yielded light curves qualitatively similar to those presented in Figure ^ and Figure 
Although one would generally expect more complicated light curves such as those depicted in 
Figure |5], a significant fraction of the geometries yield light curves in which the flux is constant 
with rotational phase when it reaches its maximum, and during the rest of the precessional period, 
the hotspot is hidden for a large portion of each rotation. 



Discussion 



We have explored several possible light curves of a freely precessing neutron star with a hotspot, 
and focussed on those whose flux is constant with rotational phase when the flux reaches a maximum 
value. These possibilities indicate that IE 161348-5055 may be a freely precessing neutron star. 
However, as the mean flux decreases from its maximum it also begins to vary with the rotation 
of the star; therefore, in the context of this model, we would expect that subsequent observations 
of IE 161348-5055 may uncover its rotational period which we would expect to be on the order of 
several seconds. The variation of the flux with the rotational phase of the star during some portion 
of the precession appears generic to freely precessing stars with a hotspot. 

The Earth undergoes a free precession with a period of about 433 days, known as the 



Chandler (1891) wobble. This is significantly longer than one would expect from the asymmetry of 



the Earth's figure due to dissipation inside the Earth (Bursa & Pec 1993). Without excitation the 
wobble would disappear within about a century, and the avenue for its excitation is still unclear. 

In neutron stars, precession has been proposed to explain long-term variations in their spin and 
pulse profiles (e.g. Davis fc Goldstein 1970| , Goldreich 197C| , Ruderman 1970| , Brecher 1972, Pines & 
Shaham 1972| , |Pines, Pcthick fc Lamb 1973 , Pines & Shaham 1974). If neutron stars rotate as rigid 
bodies the precessional period would be P/e (Pines & Shaham 1972, pines fc Shaham 1974] ) where 



P is the rotational period. |Ruderman fc Sutherland (1974)| proposed that neutron stars contain a 
super fluid component in their cores. 
I 1 

Shaham (1977) explored how the pinning of the superfluid vortices affects the free precession 
of a neutron star. He argued that the dissipation timescale for the precessional mode (t w ) is of the 
order of the postglitch relaxation time (r) times the ratio of the rotational to the precessional 
frequency, r ranges from a week for the Crab to nearly a century for 1641-45 QShapiro fc 
Tcukolsky 1983| ); therefore, depending on the nature of the superfluid coupling the precessional 
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Fig. 4. — The left panel gives the light curves as a function of precessional phase (4>) and rotational 
phase (rj) for y = 0.2 R = 11 km for a 1.4 M neutron star, k = (3 = £ = 6(1) and 7 = and 
£ = 0(1). The right panel depict the mean flux as a function of 4>. The dashed curve gives the 
best-fitting function of the form a cos <f> + b. 




Fig. 5. — The left panel gives the light curves as a function of precessional phase ((f)) and rotational 
phase (77) for y = 0. The left panel has 7 = 69°, k = 144°, j3 = 62° and £ = 115° and the right 
panel has 7 = 138°, k = 96°,/? = 126° and £ = 141°. 
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mode may last for millennia. [Shaham (1977) also found that if the star is triaxial the geometry of 
the precession is more complicated that for the purely free precession considered here. Additionally, 
the precessional frequency in this case is given by angular velocity of the superfluid component of 
the star times the fractional contribution of the superfluid to the total moment of inertia of the 
star (about one percent). 

i : — ] 

Sedrakian, Wasserman & Cordes (1999) have recently reexamined the precession of 
multicomponent neutron stars with imperfect vortex pinning and found several possibly long-lasting 
precessional modes with long periods like the precession described here. The excitation and decay 
of precessional motions in neutron stars are still uncertain. 



Evidence has been found for free precession in some radio pulsars (e.g. Cadez, Galicic fc 



Calvani 1997 , Jones 1988|) , but it is not generic (e.g. Morgan et al. 1995| ) ; therefore, the question 
arises as to which properties of a neutron star would allow or prevent it from precessing and would 
they correlate with its radio emission. Melatos (1999) and Melatos (2000) argue that precession is 



characteristic of strongly magnetized neutron stars (c.f. Equation 0). Usov & Melrose (1996) and 



Arons (1998) have proposed that strongly magnetized neutron stars are unlikely to produce radio 



emission collectively due to the formation of bound electron-positron pairs. Alternatively, Baring 



fc Harding (1997) suggest that in sufficiently strong fields (B ^ B c ) the QED process of photon 
splitting (|Adler 197l| ; |Heyl fc Hernquist 19971 ) can dominate one-photon pair production. This will 
effectively quench the pair cascade, making coherent pulsed radio emission impossible. 

Since the timescales for both the excitation and decay of precessional motion in neutron stars 
are unknown, one can appeal to the relative youth of IE 161348-5055 and the other members of 
the AXP class. They are all several thousand years old, much younger than vast majority of radio 
pulsars ( Taylor, Manchester Lyne 1993| ). The appropriate timescales for precession may simply 
be shorter than the ages of most radio pulsars while longer than those of AXPs. Furthermore, 
the hints of precession seen in the Crab pulsars ( |Cadez, Galicic Sz Calvani 1997 ) point toward this 
explanation. 

We have examined the light curves of free precessing neutron stars with a hotspot and focussed 
on those geometries which exhibit an epoch during each precessional period where the flux does 
not vary as the star rotates. These geometries account for about three percent of a random sample 
and may provide an explanation for the emission from IE 161348-5055 . If this is the case, further 
observations of the light curve from IE 161348-5055 should reveal a pulse period of the order of 
10 -4 times the precessional period of six hours. Free precession may be a hallmark of young or 
highly magnetized neutron stars, and it is a direct probe of the structure of the crust and interior 
of the neutron star and the coupling between them. 
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